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b : Abstract 



A comparative dynamical study of axial gauge QED and QCD is presented. El- 
ementary excitations associated with particular field configurations are investigated. 
Gluonic excitations analogous to linearly polarized photons are shown to acquire infi- 
nite energy. Suppression of this class of excitations in QCD results from quantization 
of the chromelectric flux and is interpreted as a dual Meissner effect, i.e. as expulsion 
from the QCD vacuum of chromo-electric fields which are constant over significant dis- 
tances. This interpretation is supported by a comparative evaluation of the interaction 
energy of static charges in the axial gauge representation of QED and QCD. 
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1 Introduction 



Confinement of colored objects is one of the fundamental properties of the strong 
interaction. Apart from the obvious importance for the spectrum of observed hadrons 
in the absence of free gluons and quarks, confinement is quite likely responsible for other 
more general properties such as the existence of Regge trajectories and the spontaneous 
breakdown of chiral symmetry. It is believed that color confinement is only one of 
the possible phases associated with the strong interaction which at finite temperature 
and possibly at finite baryon density undergoes a transition to an unconfined phase. 
Theoretical evidence based on lattice gauge calculations strongly suggests that QCD 
indeed exhibits the phenomenon of confinement and the transition to an unconfined 
phase at finite temperature. Significant experimental efforts are underway to establish 
the existence and to investigate the properties of the unconfined phase. 

Despite many efforts, a generally accepted analytical explanation or qualitative 
description of confinement within the framework of QCD is still missing. The only 
theoretical model of a confining theory, which has a well defined derivation from QCD, 
is provided by the strong coupling limit of lattice QCD. Unfortunately, this strong 
coupling limit does not distinguish QED and QCD as far as confinement is concerned. 
These two theories supposedly develop their characteristic differences only in a phase 
transition, which as a function of the coupling constant is known to occur for QED [JTJ] 
and for which no evidence is found in lattice QCD calculations. 

Canonical "gauge-fixed" formulations of QCD 0, |3| constitute a complementary 
approach to the issue of confinement. We will show how, within this framework, non- 
abelian properties of QCD provide strong signatures of confinement. Distinguishing 
elements of QED and QCD remain manifest. A variety of gauge fixed formulations 
of QCD have already been studied (cf. 0, f| ||) with the aim of clarifying the 
structure of QCD. The general strategy of all these efforts consists in explicitly resolv- 
ing the constraint equations of QCD, such as the Gauss law in the Weyl-gauge, and 
thereby eliminating the redundant variables from the QCD Hamiltonian. After this 
elimination, the Hamiltonian is, by construction, formulated in terms of gauge invari- 
ant variables (leaving aside the issue of global, residual symmetries). The resulting 
gauge fixed formulation of the dynamics is accessible to approximate treatments which 
are fully compatible with the underlying local gauge symmetry. Without gauge fixing, 
approximations generally violate gauge invariance and are therefore not appropriate 
for investigations of issues like confinement, where exact local color conservation must 
be particularly important. 

The paradigm of a gauge fixed canonical formalism is provided by the Coulomb- 
gauge representation of QED. Here, by explicitly resolving Gauss's law the longitudinal 
components of the vector potential are eliminated. The resulting description in terms 
of transverse vector fields is explicitly gauge invariant irrespective of the approxima- 
tions involved in application of the formalism to specific dynamical problems. This 
gauge choice has played a very important role also in the development of gauge fixed 
formulations of QCD 0, |, ||, |7| . However while in QED the Coulomb-gauge is singled 
out as the gauge in which static charges do not radiate, gluons and color spin remain 
coupled in QCD in any gauge. Moreover, the choice of radiation gauge variables does 
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not conform with the angular momentum algebra of the QCD Gauss law operators, 
thereby causing significant technical disadvantages. As a consequence, the Hamilto- 
nian of QCD in the Coulomb-gauge representation cannot be constructed explicitly; its 
calculation involves inversion of the Fadeev-Popov operator which can be performed 
only perturbatively. Therefore, within the canonical formalism, the Coulomb-gauge 
has been particularly useful in perturbative studies of, for example, the running cou- 
pling constant 0, |9], |1(| [11] or the small volume limit of QCD ||12|| . However dynamical 



calculations within the Coulomb-gauge are of great technical complexity . 



In the absence of compelling physics reasons, the choice of gauge has been dictated 
in most studies by formal and technical considerations. We are only at the beginning of 
a development where common properties of the resulting quite different formulations 
of QCD, such as the appearence of centrifugal barriers in the electric field energy 
13, IT5L n>|, are being recognized and their dynamical role is being understood. At 



this point it appears unlikely that the non-perturbative phenomena of QCD are most 
easily described in one unique gauge. Given the complexity of the Hamiltonian in all 
the gauge fixed formulations, it seems appropriate to attempt to identify in a first step 
those degrees of freedom whose dynamics gets particularly simplified in a special gauge 
and concentrate the studies on their dynamics. We present here the results of such an 
investigation in the framework of the axial gauge representation of QCD. 

In the axial gauge representation, the "simple" degrees of freedom are linearly 
polarized plane waves. This simplification arises by identifying one of the coordinate 
axes (the 3-axis in the following) with the spatial direction of the electric field and 
furthermore, in exploiting local gauge invariance, by identifying the color 3-axis of 



SU(2) QCD with the color orientation of the electric field. It can be shown |14| that, in 
implementing the Gauss law, all other field configurations with electric fields pointing in 
the spatial 3-direction have a longitudinal component and can therefore be eliminated. 
The simplification of the dynamics occurs for both QED and QCD. In order to control 
infrared difficulties the system is enclosed in a box with periodic boundary conditions. 

The main part of our studies will be concerned with the elementary excitations of 
these "simple" degrees of freedom. In QED, these elementary excitations are nothing 
else than photons described in this particular gauge by a single field. In QCD, the 
corresponding elementary excitations will be seen to not exhibit the standard dispersion 
equation for massless particles; rather these excitations are found to be frozen in the 
infinite volume limit. We shall interpret the disappearence of this class of excitations 
from the spectrum as an indication of the dual Meissner effect: the QCD vacuum 
apparently does not support color-electric fields which remain constant in a certain 
direction over "large" distances. The parallel treatment of QED and QCD in analogous 
gauges is very important for identifying those elements of QCD which make it differ 
so significantly from QED. Unlike lattice theories, the canonical description is not, 
by assumption, formulated in terms of compact variables; a restriction in the range 
of certain dynamical variables occurs rather as a consequence of modifications in the 
electric field energy which are typical for non-abelian theories. As another application 
of the simplified dynamics of plane wave excitations we consider the static quark- 
antiquark interaction. Again, a proper choice of coordinates and variables allows us to 
simplify significantly the dynamics and, under additional assumptions, derive within 
the canonical gauge fixed formalism the strong coupling result of lattice gauge theories. 
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2 Formalism of axial gauge QCD 



Our investigation starts with the QCD Hamiltonian in the axial gauge. For SU(2) color 
and in the presence of static charges, the Hamiltonian density is given by |I4[ 



n = tr 



E 2 , +B 2 



+ + 



Tjj^el e 3 (xj.) + ^ff 2 {x±) 



(2.1) 



In the axial gauge representation of QCD, the degrees of freedom are the (three color 
projections of the) perpendicular components of the gauge fields A± and their conjugate 
electric fields E±. (We use an unconventional sign for the electric field, so that A and 
E satisfy standard commutation relations of coordinates and momenta.) One of the 
cartesian components of the fields, the 3-component, has been eliminated. Naive gauge 
fixing gives rise to infrared problems; these can be avoided by enclosing the system 
in a box (length L) and imposing periodic boundary conditions. In this way, two- 
dimensional, i.e. ^-independent, color neutral gauge fields a 3 and their conjugate 
electric fields e 3 (x±) remain as unconstrained degrees of freedom. In turn, the two- 
dimensional neutral, longitudinal components of A± and E± are identified as dependent 
degrees of freedom and have been eliminated, i.e., 



div_i_ 



dx 3 E 3 ± 



. 



(2.2) 



(Note that the canonical commutation relations are correspondingly modified.) The 
relation between magnetic and gauge fields is the standard one, 



B 3 
B, 



d 1 A 2 -d 2 A 1 -ig [A U A 2 
d k a 3 r 3 /2 - d 3 A, 



ig 



A k ,a 3 r 3 /2 



1,2) 



(2.3) 



Except for the modifications by the two-dimensional fields, the contributions of the 
perpendicular degrees of freedom to electric and magnetic field energy are as usual; the 
second term in eq. Q2.1|) is just the contribution of the 3-component of the electric field 
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1 

1? 



dz 3 / dy 3 (1 - 8p*8 n ,o, 



^Ipq^-L? z 3 )G± pq (x±, y 3 ) 



J2wn(z 3 -y 3 )/L 



p.qjl 



^T + g(p-q)a 3 (x ± ) 



(2.4) 



which, by resolving Gauss's law, is given in terms of the other degrees of freedom, 

G ± (x) = V ± E ± (x) + ge abc ^A\ (x) (E c x (x) - ff(x ± ) 6 Cj3 ) + gp m (x) . (2.5) 



We have omitted the quark contributions in the Hamiltonian but included for later 
application static color charges as described by the density p m (x). The Hamiltonian 
density of eq. ( |2.1| ) also contains the two-dimensional, neutral color-electric field ff 



7](X±) 



d 3 y d (x ± - y±) (< 



e 3de 11 (x)E" 



[x 



Pm i X , 



which is defined in terms of the two-dimensional Green function 



1 



1 



,iq n z± 



. 2tt 

q n = -j-(ni,n 2 ) 



(2.6) 



(2.7) 
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This "electrostatic" field appears as a consequence of the elimination of the neutral, 
longitudinal components from E± and A± (cf. eq. ( |2.2| )). Finally, as a remnant of the 
Gauss law constraint, periodicity of the fields requires the neutral component of the 
charge to vanish in the space of physical states |<&), 



Q 3 |$> = 2 |d 3 x(G ± )n(x)|$) = 



{21 



Further insight into the dynamics as described by the above QCD Hamiltonian can 
be gained from a comparison with the Hamiltonian of QED in the axial gauge 



n 



1 r 

2 



E, + B 



— ° 1 9 / -» \ / — * \ 



-A 



(2.9) 



The same choice of the degrees of freedom as physical (E±) or constrained (E 3 ) ones 
preserves the similarity of the structure of the underlying Weyl-gauge Hamiltonians of 
QED and QCD. Obvious differences are connected with the color structure, e.g. in 
the definition of the QED magnetic fields (g = in eq. ( |2.3|) ), or in representing the 
3-component of the QED electric field in terms of the unconstrained variables, 



f* = 2h!o dz *!o dy ^ [l - 5nfl) 





2-nn 


2 




L 





with 



G± (x) = V ± E ± (x) + ep r ' 



[x 



(2.10) 



:2.11 s . 



In QED, the electrostatic two-dimensional field ff is generated by the density of the 
charges which again are assumed to be static, 



d 3 yd(x_L - y ± )p r ' 



[x 



[2.12) 



As in QCD, two-dimensional fields appear as a consequence of implementing the Gauss 
law on a torus. Here, the physics reason is particularly transparent. The electric 
fields e3 (x_i_) are purely transverse fields and are therefore not part of the Gauss law 
constraint. These fields together with the conjugate variables 03 (x±) describe photons 
polarized in the 3-direction and propagating in the 1-2 plane. These degrees of freedom 
cannot be eliminated while their counterparts, the purely longitudinal components of 
E±, are completely determined by the Gauss law and eliminated by the constraint 



3 Elementary excitations 

The role of the two-dimensional degrees of freedom appearing in both the QCD and 
QED axial gauge Hamiltonians can be discussed from two quite different points of 
view. On the one hand, their presence can be understood purely formally to guarantee 
consistency of the formulation and provide proper infrared behavior of the energy 
density. We note, for example, the occurrence of a 3 (x±) in the propagator defining 
£f in eq. fl2.4p . Apart from this formal role one might expect physically these special 
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two-dimensional degrees of freedom to contribute in a thermodynamic sense negligibly 
to physical observables. On the other hand, these two-dimensional fields generate 
legitimate elementary excitations of the system. In QED, any linearly polarized photon 
can, after an appropriate choice of coordinates, be described by these lower dimensional 
fields. In this sense, the degrees of freedom associated with a 3 (x±) and e 3 (x±) are of 
rather general nature, and it is only the peculiar choice of coordinates which simplifies 
significantly their dynamics. 

This simplification becomes explicit for QED by observing that after an integra- 
tion by parts in the magnetic field energy of eq. ( |2.9| ) (cf. also eq. ( |2.3| )), the two- 
dimensional degrees of freedom decouple from the remaining ones and are described 
by the Hamiltonian 

h = J d 2 x (x ± ) + | (V ± a 3 (x ± )f . (3.13) 

We compare this Hamiltonian with the one for the corresponding elementary exci- 
tations in QCD. Here, the color dynamics couple the two-dimensional fields to the 
perpendicular degrees of freedom. In QED, the presence of charged matter would also 
induce such a coupling. In a first step we shall neglect this coupling and keep - with 
the electric field energy - only the abelian contribution to the magnetic field energy. 
The resulting Hamiltonian for a 3 (xj_) of QCD is, by construction, identical in structure 
with the QED one of eq. ( |3.13|) but for the missing hermiticity of the QCD electric 
field operator e 3 (x±) (cf. eq. (|2.1|)). This difference has important consequences. The 
hermiticity defect arises since SU(2) group elements (W), actually loops around the 
torus along the 3-direction, have been parametrized in terms of elements of the algebra 
(as), 

W = exp {igLa 3 f/2} = cos (gLa 3 /2) + id 3 f sin (gLa 3 /2) . (3-14) 

This is analogous to the hermiticity defect of the radial momentum operator when 
using polar coordinates and also leads here to modifications in the kinetic energy. In 
the Schrodinger representation, expressing the electric fields in terms of functional 
derivatives, one finds 

1 5 5 

e\ (zj_) e 3 (x ± ) = — - — T—r- — — -J (a 3 (x±)) - — — - . (3.15) 
J (a 3 (x ± )) 5a 3 [x±) 6a 3 (x ± ) 

The Jacobian J (a 3 ) is the Haar measure of SU(2), 

.2/1 



J (a 3 (x±)) = sin y-gLa 3 (x±)j (3.16) 

and thus also appears in the volume element when evaluating matrix elements in terms 
of wave functionals of a 3 (x±). As the transformation to polar coordinates becomes 
singular when attempting to define a direction for a vector of vanishing length, singu- 
larities in the kinetic energy eq. (|3.15|) arise for gLa 3 = 2im in the parametrization 



of eq. ( |3.14|) . It is convenient to transform the kinetic energy into the standard form 
by defining, in analogy to the Schrodinger equation in polar coordinates, appropriate 
"radial" wave functions. We introduce for this purpose a lattice in the 1-2 plane (lattice 
constant £) and define the (rescaled) variables at the lattice sites b = x±/£, 

n =lgLa 3 (b£) . (3.17) 
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The radial wave functions $ are defined in terms of the original wave functions, which 
are the projections of the physical states |<&) (cf. eq. ( |2.8|) ) onto the field eigenvectors 
\(p), in the standard way 

<%] = <%]nisin(^)| . (3.18) 
b 

In this way formally identical "free" Hamiltonians of the two-dimensional gauge degrees 
of freedom (p for QED and QCD are obtained 

h = h e + h m , (3.19) 

with the electric and magnetic field energies given by 

h* = _^V^ 
8& ^ dipl ' 

b Y b 

h m = -|fE {Vb + 6-nf- (3-20) 

y b,8 

The fundamental vectors of the lattice (1,0) and (0,1) are denoted by 5. Although 
of the same structure, the Hamiltonian ( |3.19p acts on wave functions belonging to 
different spaces in QED and QCD. The space of wave functions is restricted in QCD 
by the constraint 

=0 whenever cp? = nn for some b . (3-21) 

In the transformation of the Hamiltonian to "radial" variables no effective potential 
but an energy shift 

E^E-££ (3.22) 
appears, which has been suppressed in eq. (|3.20|) . 



3.1 Photons in the 1-2 plane 

In order to display the non-trivial dynamics described by the Hamiltonian h in conjunc- 
tion with the constraint of eq. (|3.21|) we consider first the case of electrodynamics. As 



is well known, in this case eigenstates and energies are determined by discrete Fourier 
transformation 

n=^H^ rk/K n- (3-23) 

k 

The number of degrees of freedom is K 2 with 

K=j (3.24) 

and limits the sum in eq. (|3.23| ) to \ki z\ < K — 1. In this way, the dispersion relation 
for the (lattice) photons, 

4 „ . „ / ixdk 



UJ 



2 _ ^ „;„2 
k 



£ sin 2 (^] (3.25) 
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is obtained, with the continuum limit 

2 



2-1^1 for \k\ « K . (3.26) 
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Thus the two-dimensional field (x±) describes free photons propagating in the 1-2 
plane with polarization in the 3-direction. We also note that the dependence of the 
Hamiltonian on the coupling constant is of no relevance, a simple rescaling 

W^gLWtpj; (3.27) 

eliminates all but the dependence on the lattice spacing I. 

In momentum space, the ground state wave function factorizes into the contribu- 
tions of the normal modes 



m = n 

k 



U 2 uA l/2 ( M 2 _ _ N 
jf) exp —^T u k<Pk<P-k 



(3.28) 



ng 2 L J V g 2 L 

and thus describes highly correlated zero-point motions in configuration space. 

3.2 Jacobian and structure of the Hilbert space 

The nature of the elementary excitations as described by the Hamiltonian ( |3.19| ) is 
qualitatively altered by the presence of the constraint (|3.21|) on the wave function. 
This difference in dynamics is easily appreciated in a mechanical interpretation. While 
the Hamiltonian (|3.19|) , ( |3.20|) describes a two-dimensional system of mass points which 
are coupled harmonically to the respective nearest neighbors with the photons as the 
normal modes, the constraint ( |3.21| ) can be visualized as walls of an infinite square 
well which limit strictly the amplitude of oscillations of the individual mass points. 
Apparently, the character of the normal modes of this system depends on the relative 
importance of the harmonic nearest neighbor couplings and the constraining force. The 
relevant parameter which controls the dynamics is the ratio of constants multiplying 
electric and magnetic field energy, 

k = g 2 ^ . (3.29) 

Unlike the case of electrodynamics, this ratio receives significance by the constraint on 
the wave function which prevents the kinetic energy from becoming arbitrarily small 
for certain field configurations. 



The presence of the constraint ( 3.21|) on the wave function not only changes signif- 



icantly the physical properties of the system but also requires, in comparison with the 
abelian case, a very different method of solution. As a consequence of the constraint, 
the probability current between regions of the wave functions separated by zeroes of 
the Jacobian vanishes. Therefore the system defined by eqs. ( [3.19| ) - ( |3.21f ) possesses 
an infinity of conserved charges 

Ql = / d<pS (y - t ) (3.30) 
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(here we have made explicit the operator character of ipjf). These charges specify the 
probability of the system to live in one of these intervals. In general, stationary states 
do not exist for an arbitrary choice of these charges; they rather require each of these 
charges to be non-vanishing only in one of the intervals. For illustration we consider 
the case of two lattice sites. 



In accordance with the constraint ( 3.21|) we decompose the wave function 



1>(<Pl, <&) = E &n 1 M®n 2 (<P2)lpn 1 ,n 2 (<Pl,<P2) , (3.31] 



ni ,n2=— oo 

where G n denotes the projection on the fundamental intervals 

Most important, the wave function constraint decouples the solution of the Schrodinger 
equation in the different intervals; consequently eigenvalues and eigenf unctions can be 
obtained by choosing the wave function to be non- vanishing in one of the intervals only, 
i.e., 

^m^^b^^O only if m = n°, n 2 = n° . (3.33) 

The complete spectrum is determined by variation of n\ and n®. In general, the po- 
tential energy and therefore the spectrum depends on the choice of the intervals and 
therefore the system cannot coexist in different intervals. When symmetries are present 
however, the wave function corresponding to identical eigenvalues in different intervals 
can be arbitrarily distributed over these intervals. For example, this happens as a 
consequence of the displacement symmetry of the Hamiltonian ( |3.19| ), ( |3.2(J| ); a global 
shift of the intervals 

n i ~* n l + m ( 3 - 34 ) 
does not change the spectrum and the wave functions can be linearly superimposed. 
Further possibilities for combining wave functions arise, for instance, by the following 
more symmetric form of the discretized magnetic field energy, 

hm - ^E (1 " 008 2(^-^5)) . (3.35) 

3 b,S 

The presence of symmetries makes the definition of the eigenstates ambiguous; however 
the theory contains no operators which would connect different intervals with each other 
and thereby be sensitive to these ambiguities. In the following we shall restrict our 
calculations to one definite (site-independent) choice of the fundamental intervals, 

Ql = 5n,o ■ (3.36) 

The site-independence guarantees that the state of lowest energy occurs in this sector 
of the Hilbert space (with the parametrization ( |3.35| ) of the discretized magnetic field 
energy, the spectrum is independent of the choice of the fundamental intervals). 
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3.3 Gluonic excitations in the 1-2 plane 



The strong coupling limit (k ^> 1) is physically most important. We note that this limit 
is determined by the strength of the dimensionless coupling constant g relative to a 
negative power (—1/6) of the number of degrees of freedom in the system (cf. (|3.29| )). In 
the absence of a dimensionful quantity, the continuum (£ — » 0) and "thermodynamic" 
(L — > oo) limits are not distinguished. We therefore expect the dependence of the 
coupling constant g on the number of degrees of freedom to be dictated by asymptotic 
freedom, 

S 2 (L/£)cxl/m(L/£) , (3.37) 

and the thermodynamic limit to correspond to the strong coupling limit. In this limit, 
the electric field energy dominates and the stationary states are simply given by exci- 
tations of the degrees of freedom at the individual lattice sites. The wave function of 
a stationary state can be written as 



*wM =n 

b 

with the energy eigenvalue 



sin [npp } 



(3.38) 



£« = §EH|- ( 3 - 39 ) 

b 

In particular, in the ground state all the individual degrees of freedom are in the lowest 
n = 1 state, and the value 

n 2 T 3 

Eo = ?-=r (3-40) 

is obtained for the ground state energy. This is just the negative of the energy shift 
associated with the transformation to radial coordinates (cf. eq. ( |3.22|) ). Thus the 
total ground state energy is actually zero and the corresponding ground state wave 
function is an eigenfunction of the modulus of the chromo-electric field operator with 
vanishing eigenvalue (up to singular contributions arising from possible discontinuities 
at (p — nn). States of lowest excitation energy are obtained by exciting a degree of 
freedom at one particular site into its first excited state, with excitation energy 

Corrections to the strong coupling limit can be calculated perturbatively. First 
order perturbation theory in h m (cf. eq. (|3.20|) ) yields for the ground state energy 

^ g 2 L 3 AL fir 2 \ . . 

E » = -w + w{^- 1 ) ' (3 ' 42) 

The result justifies the perturbative treatment of the magnetic field energy for suffi- 
ciently large number of degrees of freedom k 2 = g 4 (L/A£) 2 ^> 1. 

Excited states of the system are highly degenerate in the strong coupling limit. This 
degeneracy is lifted due to the magnetic coupling of the degrees of freedom at different 
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lattice sites. We consider the energetically lowest excitation at the site b described by 
the wave function 



2N l/2 
7T 



sin ( 2 ^) n 



2 x 1/2 



7T 



sin <y9p 



(3.43) 



In the subspace of these excitations, the nearest neighbor magnetic coupling becomes 
diagonal for the eigenstates of the translation operator (6 — > b + 5) 



b 

The calculation of the excitation energies is straightforward; we find 



(3.44) 



AE 



3 g 2 L 64 1 



k 8 I 2 9n 2 g 2 L 



sin 



(3.45) 



corresponding to the state Q3.44Q labeled by the momentum k. 

The above analysis of the strong coupling limit is the central part of our investiga- 
tions. We now present the physics implications of the results. As already emphasized, 
the ground state of the system is, in the strong coupling limit, an eigenstate of the 
electric field operator with vanishing eigenvalue. This possibility for a ground state 
with vanishing ^-independent electric fields arises due to the Jacobian in the kinetic 
energy. In QED such states are not normalizable and would entail infinitely large fluc- 
tuations in the magnetic field energy. Thus the structure of the vacuum concerning 
these ^-independent fields is very different in the abelian and non abelian theory. The 
virial theorem applies to QED: magnetic and electric fields contribute equally to each 
normal mode. Hence, the expectation value (E 2 — B 2 ) vanishes. In QCD the Jacobian 
invalidates equipartition. Chromoelectric ^-independent fields are absent; in turn, the 
fluctuations in the magnetic field at different lattice sites are not correlated and the 
ground state energy is due exclusively to these uncorrelated magnetic field fluctuations. 
Similarly, the "gluon condensate" is dominated by the magnetic field contribution, 



(E 2 -B 2 ) = -(B 2 ) <0 . 



(3.46) 



It is interesting that even in the crude approximation of keeping only one particular 
kind of gluonic degrees of freedom, this model displays features which are reminiscent 
of the phenomenology of the "magnetic QCD vacuum" . 

Concomitant with these qualitative differences in the structure of the ground states 
is the very different nature of the elementary excitations. Built on the highly corre- 
lated QED ground state the photons appear as collective excitations with excitation 
energies vanishing in the long-wavelength limit. In QCD, the elementary excitations 
are localized in configuration space and are due to formation of non-vanishing electric 
flux. Intuitively, the expression ( |3.41| ) for the excitation energy can be seen as a result 
of a quantization of the chromoelectric flux. In the absence of couplings to the other 
degrees of freedom, the chromoelectric fields formed in the elementary excitations with 
lowest excitation energy are located on just one transverse lattice site. Thus the flux 
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tubes are infinitely thin (i.e., the area is oc £ 2 ) and extend over the whole system in 
the 3-direction. The chromoelectric flux $g associated with each component of the 
standing waves (eqs. (|3.43|) , (|3.44j) ) is quantized, 



$ B = ng . (3.47) 

(Note that since the standing waves are not eigenstates of the electric field operator, 
one has to decompose them into travelling wave components in order to "see" this flux 
quantization.) Beyond the strong coupling limit, waves propagate through the medium 
and transport the electric flux across a transverse plane. The spectrum acquires a band 
structure with a band width oc l/(g 2 L) (cf. eq. ( |3.45| )). Threshold energy and mass 
of the associated "particles" are cx g 2 L/£ 2 and tend to infinity with increasing size 
of the system. Irrespective of the value the transverse momentum, these elementary 
excitations never approach the free gluon limit. 



4 Gluonic couplings 

Here we continue our investigation of the elementary excitations of QCD as described 
by the two-dimensional fields (cf. eq. ( |3.17| )). Unlike in QED, where these degrees 
of freedom describe non-interacting photons in the absence of charged matter, in QCD 
such a decoupling does not occur. Here we study - within perturbation theory - the 
effect of the other degrees of freedom of QCD on the dynamics of the yg. We shall show 
that the effect of the infinity of other degrees of freedom is not sufficient to overcome 
the dominant role of the kinetic energy of these particular two-dimensional fields. The 
starting point of our investigation is the Hamiltonian of eq. ( |2.1| ) applied to SU(2) 
QCD. With the coupling constant g treated as a small parameter, the neutral gluons 
as described by the conjugate pair of vector fields A 3 ± , E\ decouple from the other 
degrees of freedom and will not be considered further. At this point it is convenient to 
define "charged" gluon fields 




(4.48) 



The Hamiltonian to be considered in the following, 

h = h e + h m + h'[<p] , (4.49) 

contains apart from the already discussed electric and abelian magnetic field energy 
the Hamiltonian h'[ip] of the charged gluons coupled to the two-dimensional degrees of 
freedom ip^, 

h'[ip] = Jd 3 x(nl{x)n ± {x)+f3\x)(3{x)) (4.50) 

L f f L f L - -a. i2nn(z-z')/L ^ ^ 

+ - U 2 x ± dz dz'Y,v±n±{z±,z)-, 7 zrvaV±n±@±,J) ■ 

4 J JO JO n [TTTl — CpyXx)) 

The kinetic energy of the charged gluons receives contributions from the first and last 
term of the r.h.s. of the above equation. The last term also acts as a centrifugal 
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barrier on the two-dimensional degrees of freedom (p(x±) (here, we have adopted a 
continuum notation for these variables, since most of the calculations can be carried 
through without resorting to the lattice in the perpendicular direction). The magnetic 
field (3 is defined as 

p 3 = d x ® 2 - d 2 ® x , p 2>1 = ± (da - 2icp(x ± )/L) $ li2 . (4.51) 

As the centrifugal term, the magnetic field energy contains a coupling between charged 
and neutral two-dimensional gluons via the covariant derivative (83 — 2i<p/L). 

The singular coupling (cf. eq. ( f4.50|) ) between the two-dimensional fields <p(x±) 
and the charged gluons prevents straightforward application of perturbation theory. 
Expansion of the centrifugal term in eq. ( |4.50| ) in terms of the variables <f{x±) yields 
non-integrable, infrared singularities. Thus the following considerations are more gen- 
erally required for a systematic expansion in g when these two-dimensional fields with 
their explicit non-perturbative dynamics are involved. 

To define an appropriately decoupled zeroth order Hamiltonian we introduce an 
x_i_-independent external, neutral field \ which represents the average equilibrium po- 
sition of the ip(x±). In this way a well defined, though parameter dependent interaction 
Hamiltonian is introduced. Eventually, the external field X w iU be treated as a varia- 
tional quantity, used to minimize the ground state energy. Formally, we write 

h'[ v } = h'[ X ] + (hy]-h'[ X ]) (4.52) 

and define the ^-potential energy resulting from the coupling to the charged gluons as 

h b *=(& x \(h'[p]-h , M)\Q° x ) . (4.53) 

The state vector |<&°) denotes the x-dependent ground state of the non-interacting 
charged gluons, 

h'[x]K) = E omK) ■ (4-54) 

In the first step towards the determination of h- m t, we diagonalize the free charged 
gluon Hamiltonian. Since h'[x] is quadratic in the charged gluon field operators, this 
diagonalization is achieved by a Bogoliubov transformation and yields the result 

h'bc] = E IK") I + a[(-n)a x (-n) + l) . (4.55) 

The operators a\(n), P\(n) create the two (SU(2)) charged massless gluon states with 
polarization specified by A and momentum given by k(n), 

2n 

^(n) = — (ni,n 2 ,n 3 ) , (4.56) 
while their energies are determined by the vector 

p{n) = — (ni,n 2 ,n 3 - xM ( 4 - 57 ) 

associated with the covariant derivative. The coupling to the external, spatially con- 
stant neutral gluon field X induces an asymmetry in the energies of the two charge 
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states. Furthermore, the independence of the energies \p(n)\ of the polarization A re- 
quires a consistent treatment of both magnetic and electric coupling of x to the charged 
gluons. For vanishing \ and with the n = term in the centrifugal term (cf. eq. (|4.50|) ) 
disregarded, h'[x] represents (apart from the corresponding two-dimensional fields) the 
free Hamiltonian for two types of photons in the axial gauge representation. By the 
Bogoliubov transformation this Hamiltonian gets transformed into the Coulomb-gauge 
representation and thereby the non-locality in the electric field energy is eliminated 
T7|. Having determined the Bogoliubov transformed ground state it is straightforward 



to calculate the separate contributions to the ground state energy. The color-electric 
contribution hf nt to the potential energy h mt (cf. eq. ( |4.53| )) arises from the centrifugal 



barrier term in eq. ( f4.50| ), the color-magnetic one h™ t from the perpendicular compo- 
nents of the magnetic field, 

h int = ht nt + hZ . (4.58) 

To evaluate these matrix-elements we replace at this point the corresponding integrals 
over the perpendicular coordinates by sums and obtain 



"■nit 




^/L\ 2 A kl 



Pn3 + kn± 



* = ^sHCTfJ-^ftrf • (4 - 59) 



b,ft 

(Here, we have used the shorthand notation k n 3 = k^n), k n ±_ = k±(n), etc.) As could 
be expected, these energies are divergent despite the subtraction of the most singular x 
and (^independent terms in the definition ( |4.53| ) of h int . Using a heat kernel regulator 
of the form 

e-VPn3+ fc L (4.60) 
and identifying the inverse momentum cut-off with the lattice spacing 

A = 7 £ , (4.61) 

we obtain to leading order in £/L 

= ^E{2(x-^)oat^+(x- W ) 2 ^} > 
/C = ^E{2(^-x)cot X+ (^- X ) 2 -i-} . (4.62) 

Determination of the value of the external variable x requires computation of the 
expectation value of h int in the unperturbed ground state (rag = 1 in eq. ( |3.38| )). It is 
easily seen that the interaction energy becomes minimal at 

X = vr/2 , (4.63) 

i.e. at the midpoint of the intervals in which the (p% degrees of freedom move. This 
choice also minimizes the total energy; the zero point energy in h'[x] (eq. ( |4.55p ), 
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though x-dependent, does not contribute to leading order in IjL. Thus our final result 
for hi nt reads 

hi* = ^ E { (tt/2 - n ) 2 (l + (sin ^) " 2 ) + 2 (tt/2 - n ) cot ^} . (4.64) 

The properties of the potential energy h int are directly correlated with the structure 
of the centrifugal term and the non abelian magnetic field energy in eq. ( 4.50Q . This 



potential energy is local, i.e., it does not connect degrees of freedom at different lattice 
sites with each other. It therefore does not significantly change the formation of waves 
across the transverse plane nor their suppression in the strong coupling limit (« ^ 1). 
With the centrifugal term, h int is singular at (p? = 0, it. Close to the origin the kinetic 
energy of the "radial" motion (cf. eq. ( p.20| )) is thus supplemented by a centrifugal 
barrier 

d 2 d 2 n£ 1 

-^p-^! + 7?z^r (4 - 65) 

The presence of the centrifugal barrier forces the full wavefunction ($ in eq. ( 3.18Q ) to 



vanish at the origin. In addition to the centrifugal barrier, the potential energy ([4.64 ) 



exhibits a repulsive Coulomb-like singularity. The modification of the ground state 
energy due to h int is given by 

Thus coupling to the charged gluons gives rise to a non- vanishing and, in the continuum 
limit, divergent energy density (energy per unit area). This is the leading contribution 
in the continuum limit (taking into account the redefinition of the energy in eq. ( |3.22| )). 
Nevertheless, the dominant operator of the Hamiltonian ( 4.49|) remains the electric field 
energy h e (cf. eq. ( [3.201 )). 

Finally we discuss the result ( |4.63| ) concerning the average field which, in view of 
the singularities at 0, n, is plausible. This result actually reflects a symmetry of the 
exact Hamiltonian which guarantees quite generally that the potential energy of the <£>g 
degrees of freedom must be stationary at tt/2. As shown in [[14]] , the Hamiltonian (|2.1|) 
is invariant under "displacements", "central conjugations", and reversal of the color 
3-axis. Following arguments developed in the context of 1+1 dimensional QCD JTB| 
we can combine these transformations to a relevant symmetry transformation which 
changes the variables in the following way 

<S ± (x) -> e - 2i7rX3/i $i (x) , n_L (£) -> e 2lnx " /L li[ (x) 

Al (x) -> -Al (x) , El (x) -> -El (x) . (4.67) 

The invariance of the exact (eq. ( |2.1| )) as well as of the approximate (eq. ( |4.49| )) 
Hamiltonian under this transformation is easily verified. After integrating out the 
charged gluons, the only remnant of the above symmetry is the transformation of the 
variables <^g- The symmetry therefore reduces to the invariance of h- int under a common 
reflection of the variables y?g at 7r/2. This choice in turn implies a charge conjugation 
symmetric Hamiltonian (cf. eqs. (|4.55| ), ( |4.57| )). As the net effect of the coupling to 
the two-dimensional fields, the spatially periodic charged gluon fields acquire energies 
\p(n) | which actually correspond to anti-periodic boundary conditions. 
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5 Physics in axial gauge QCD 



On the basis of the results derived in the preceding sections, we shall develop in this 
concluding section the physics picture of QCD as it presents itself in the axial gauge. 
The essential result of our investigations is the suppression of plane wave type exci- 
tations, i.e., excitations where the electric field is constant along a particular spatial 
direction. In our choice of coordinates and gauge, we have identified the spatial 3-axis 
with the direction in which the field does not vary and have used the local gauge invari- 
ance to have the field pointed in the color 3-direction everywhere in space. With this 
choice of coordinates, these particular degrees of freedom - neutral, linearly polarized 
gluons - are described by two-dimensional fields resulting in a significantly simplified 
description of the dynamics. Such a simplification is not too surprising considering the 
fact that in QED any linearly polarized photon, with the corresponding choice of coor- 
dinates and in the absence of matter, is described by a non-interacting two-dimensional 
scalar field. In QCD the decoupling from the other degrees of freedom is not complete, 
as happens in QED when coupled to charged matter. A characteristic difference to 
QED is that the gluon self interactions give rise to a significant modification of the 
kinetic energy (i.e. electric field energy) of these particular degrees of freedom. This 
modification is a direct consequence of the projection onto physical states as obtained 
in the implementation of the Gauss law. It is analogous to the modification of the ki- 
netic energy of a quantum mechanical particle when projected onto states with definite 
angular momentum. As a result of this change, eigenstates to this modified kinetic en- 
ergy operator of QCD become normalizable, i.e., the electric flux associated with these 
normalizable states gets quantized. Excitation energies of states with non-vanishing 
electric flux diverge linearly with the size of the system. Up to this point, the size 
of the system has been treated as a formal parameter introduced to properly define 
the theory in the infrared. Whenever L is large compared with the relevant intrinsic 
length scale - 1 /A after renormalization or, phenomenologically, the size of hadrons - 
finite size effects should be negligible. Success of present lattice calculations actually 
indicates that the size of the system may not have to be much larger than hadronic 
sizes. In this sense, our results suggest that the QCD vacuum does not support color- 
electric fields which are constant over distances larger than hadronic sizes. Excitation 
energies associated with such fields grow with the distance over which such fields are 
constant. This behavior is independent of their variation in transverse directions, i.e., 
independent of the associated total momentum. Thus finite excitation energies can re- 
sult only for structures radically different from plane waves. In summary, color-electric 
fields constant over large distances are expelled from the QCD vacuum, i.e. the QCD 
vacuum exhibits the dual Meissner effect. This in turn implies that magnetic fields 
cannot be correlated at points in space separated by distances over which constant 
electric fields cannot extend. 

As an application of this picture, we consider the interaction between two static 
color charges in comparison with the electrostatic interaction in the Maxwell theory. In 
the axial gauge representation of QED, the Coulomb-interaction between static charges 
is not manifest, rather a linearly rising potential appears in the expression ( 2.10 ) for 



the field energy. This contribution is however only part of the electrostatic interaction. 
In the axial gauge representation of QED (as in any but the Coulomb-gauge) static 



15 



charges couple to the radiation field as described by the interference between V±E±_ 
and p m in eq. (|2.10 ). It is easy to decouple the radiation field from the charges by 



shifting the electric field 

E ± ^E ± + eV ± ^p m . (5.68) 

Thus by introducing a pair of static charges (with net charge zero at a distance d 
from each other) into the QED vacuum, the equilibrium position of the electric field 
changes. The radiation field oscillates around the electrostatic field generated by the 
static charges. For harmonic oscillations, such a shift in equilibrium position is irrel- 
evant apart from a "c-number" change in the zero point energy, which is nothing else 
than the electrostatic Coulomb-energy 

u = ss ■ < 5 - 69 > 

In QCD a corresponding shift can also be performed; it decouples however "static" 
color charges and radiation field only perturbatively. Obviously, the dynamics of the 
color spin of the static quarks necessarily remains coupled to the gluons. On the other 
hand our above discussion suggests that the QCD vacuum resists such a shift in the 
color-electric radiation field. We have seen that the vacuum does not support electric 
fields which have small variations over large distances. Such fields are introduced 
however in the shift ( 5.68|) . In order to follow more closely how a linearly dependent 



potential between static quark and antiquark arises in axial gauge QCD we assume - 
beyond our above considerations - that the QCD vacuum is void of any color-electric 
field. Introducing a static quark- antiquark pair requires adjustment of the vacuum 
to account for the correspondingly modified Gauss law. In a gauge fixed formalism, 
this readjustment is enforced by the "color-electrostatic" field energy appearing in the 
Hamiltonian in the process of implementing the Gauss law. In the axial gauge, the 
energy density of the electrostatic field is given by the operator | (£| + fj 2 ) in the 
Hamiltonian of eq. (|2.1|). With an appropriate choice of the coordinates, it is possible 
to account for the necessary readjustment of the vacuum by changes in one of the 
degrees of freedom only. For this purpose, the coordinates have to be chosen such 
that quark and antiquark are located on the 3-axis. In this case, we may assume that 
the color-electric fields E±(x) remain undisturbed and essentially zero. Also, due to 
the neutrality condition which the static charges have to satisfy, the two-dimensional 
neutral, color-electric fields ff(x±) remain zero. Thus introduction of a static quark- 
antiquark system affects only the neutral two-dimensional gluons ip^ in the Hamiltonian 
(|2.1| ). Furthermore, as we have seen, degrees of freedom at different lattice sites are 
decoupled in the strong coupling limit. Only the neutral degree of freedom p^ at the 

site 60, corresponding to the transverse coordinate of the quark- antiquark position, is 
coupled to the static charges. The Hamiltonian describing this system of quantum 
mechanical degrees of freedom reads 



Sh= -ww+ u - (5 ' 70) 



60 
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The coupling potential is the electrostatic field energy £f/2 (cf. eq. Q2.4Q ), 



u 



g 2 £ 2 L fL 



f dz 3 f dy 3 Y, (1 - $p.q$n,o) 



p™Mb ,z 3 )p™Jb Q ,y 3 ) 



D i2wn(z 3 -y 3 )/L 



nn+ (p-q)<p(b Q ) 



(5.71; 



which is given in terms of the static color charge density 



P. 



pq 



(M) 



9 

2£ 2 



a=q,q 



i>f (b, z) r P if, z) - -f £ c f z) r r (f 

1 r=l,2 



(5.72) 



Here we have represented the color degrees of freedom of the static quarks in second 
quantized form. The operator ip^ (b, zj creates a quark (a = q) or antiquark (a = q) 

at the position ((b£, z)). With our sequence of approximations we have reduced the 
calculation of the static quark interaction energy to a problem of coupled quantum 
mechanical degrees of freedom. The above Hamiltonian 5h is, after redefinition of the 
coupling constant, identical with the Hamiltonian describing the interaction of static 
charges in SU(2) QCD in one space dimension. The spectrum of this Hamiltonian has 
been determined in |TJ5|] and yields for the interaction energy of static color charges 



U 



3g 2 d 
Si 2 



(5.73) 



This result agrees with the strong coupling limit of the SU(2) static quark-antiquark 
interaction in Hamiltonian lattice gauge theory | |20|| . It is quite satisfactory that ap- 
proximations with similar physics content, the common assumption of vanishing color- 
electric fields, yield in the two quite different formal approaches the same physical 
results. This is particularly remarkable in view of the different ways in which the 
gauge symmetry is treated. 

For our derivation, the proper choice of coordinates has been instrumental. It is 
instructive to consider the problem with a choice of the 3-axis which does not coincide 
with the direction of the quark-antiquark dipole moment. In this case, the color charges 
of quark and antiquark contribute to the y 3 , z 3 integrations for different values of the 
perpendicular coordinates b. As a consequence, in the corresponding n = terms, 
quark and antiquark contributions do not cancel. They rather give rise to a centrifugal 
barrier with strength oc g 2 L/£ 2 . In the continuum limit, an infinite interaction energy 
results. Thus, with this choice of coordinates, the assumption of vanishing color-electric 
fields E± is unattainable. The system has to respond to the non-vanishing net charges 
arising in the x 3 integrations and compensate these charges by color-electric fields. 
Thereby a flux tube between quark and antiquark is formed which, with such a choice 
of coordinates, involves the complicated dynamics of the coupled E± fields rather than 
the simple dynamics of the single degree of freedom (p^. 

Comparison of the result ( |5.73| ) with the expression ( |3.41| ) for the energy associated 
with the lowest excitation of the two-dimensional degrees of freedom ip? supports the 
arguments concerning the quantization of electric flux (cf. eq. ( |3.47D ). We observe 
that unlike the lattice formulation, the axial gauge representation does not introduce 
compact variables which by construction lead to such a quantization of the electric flux. 
In axial gauge QED no such phenomenon occurs. It is the appearence of the Jacobian 
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which forces the relevant QCD degrees of freedom not to be periodic but rather to 
be constrained to a compact interval which is defined by consecutive zeroes of the 
Jacobian. Thus most crucial properties of the two-dimensional degrees of freedom 
pertinent to the issues of confinement and the dual Meissner effect are traced back 
to a fundamental difference between QED and QCD which becomes manifest in the 
gauge fixed formulation. In the lattice formulation of QCD, a similar dynamical role 
is played by the Haar measure in the definition of the partition function. Our findings 
are reminiscent of the disappearence of confinement when the Haar measure is replaced 



by a constant [21 



In our discussion of the interaction energy of static quarks, a series of approxima- 
tions has been necessary whose validity might be difficult to assess quantitatively. On 
the other hand, our treatment of the dynamics of the two-dimensional degrees of free- 
dom ip has been quite straightforward and therefore provides a novel access to the 
strong coupling limit of QCD. As the fundamental approximation in this approach, the 
coupling of the neutral <p (x±) to the charged gluon fields has been neglected. With this 
coupling taken into account, the flux tubes are expected to acquire finite transverse 
extension and thereby the strong coupling result for the excitation energy (|3.41 ) to be 
determined by the string tension 

3 -f - (5-74) 

We have performed a first step beyond the "strong coupling limit." The singular nature 
of the coupling to the charged gluons prevents straightforward application of pertur- 
bation theory and we have determined in a variational calculation the average effect of 
the charged gluons on the (p (x±) ground state. This improved ground state can serve 
as starting point of a systematic expansion. Calculation of this lowest order correction 
reveals the presence of two different formal, small parameters involved in such an ex- 
pansion. On the one hand, ordinary perturbative treatment of the non-linearities of 
QCD in both the non-abelian part of the magnetic field energy (cf. eq. fl2.3| )) as well 
as in the non-abelian contribution to the perpendicular Gauss law operator G± ( |2.5|) 
naturally requires 

g < 1 . (5.75) 



On the other hand, as the comparison of the correction ( [4. 66 ) with the lowest order 



ground state result (|3.40p shows, such an expansion involves as a small parameter 

JL « 1 ' < 5 ' 76 > 

Obviously these two expansions can be made compatible only if the coupling constant 
g is considered a function of the number of degrees of freedom (L/l). With the depen- 
dence 

g 2 (x l/ln(L/£) (5.77) 



expected to hold in the continuum limit, the two requirements (|5.75 ), ( 5.76|) are indeed 



compatible. It is encouraging that the non-perturbative effects indicating confinement 
require a behaviour of the coupling constant in the thermodynamic limit which is 
compatible with asymptotic freedom. It is tempting to speculate that absence of con- 
finement in Higgs-gauge theories (e.g. SU(2) with a scalar doublet) is in turn related to 
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a behaviour of g in the continuum limit of such models [22| which, in contradistinction 
to Q5.77 ), is controlled by a cutoff rather than by the number of degrees of freedom. 
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